Infrared behavior of the Faddeev-Popov operator in Coulomb gauge QCD 
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We calculate the eigenvalue distribution of the Faddeev-Popov operator in Coulomb gauge QCD 
using quenched SU(3) lattice simulation. In the confinement phase, the density of the low-lying 
eigenvalues increases with lattice volume, and the confinement criterion is satisfied. Moreover, even 
in the deconfinement phase, the behavior of the FP eigenvalue density is qualitatively the same as 
in the confinement phase. This is consistent with the fact that the color-Coulomb potential is not 
screened in the deconfined phase. 
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C\J \ I. INTRODUCTION 

x> : 

One of the most challenging issues in elementary particle and nuclear physics is to understand the confinement of 
the quarks and gluons in quantum chromodynamics (QCD). To understand the mechanism of confinement, there are 
j several approaches in which topological objects are responsible for color confinement. Color monopoles in maximal 
Abelian gauge and center vortices in maximal center gauge are well-known examples [J • To clarify the mechanism of 
confinement, it is important to choose a proper gauge to extract the relevant degrees of freedom for color confinement. 

The confinement mechanism in Coulomb gauge QCD has received a lot of attention recently. Coulomb gauge is a 
physical gauge in the sense that there are no unphysical degrees of freedom. Accordingly, Coulomb gauge is convenient 
for a variational approach and many attempts have been made to examine color confinement and hadron spectroscopy 
fSj ' [3, BL H, IE 0JH- It is firstly discussed by Gribov in '70s that the instantaneous interaction provides the long-range 
interaction [3] , and this is further elaborated by Zwanziger recently fic| . In Coulomb gauge, the time-time component 
f*"*- ■ of the gluon propagator can be decomposed into the instantaneous part and the noninstantaneous (polarization) part 

O; a 

D 00 (x,t)=I(x)6(t)+P{x,t). (1) 

Oh, It has been found that Dqq is a renormalization-group invariant and, as a result, both / and P are separately 
JJ^ ■ renormalization-group invariants [l(| HH, 03 • Renormalizability of Coulomb gauge QCD was also studied within the 
Hamiltonian formalism [lol [T3 | and the Lagrangian formalism [l3j . Furthermore, Zwanziger showed that the color- 
Coulomb potential which is the instantaneous interaction energy between heavy quarks is stronger than a physical 
' potential: 

S : Vcoui(#) > V(R). (2) 

_ _ i 

This inequality tells us that the necessary condition for the physical potential being a confining potential is that the 
color-Coulomb potential is also a confining potential, i.e., "no confinement without color-Coulomb confinement" [lil ]. 

In SU(2) lattice calculations, it was reported that the instantaneous part of the gluon propagator, I(k), is strongly 
enhanced at k = [SEH- Furthermore, the recent Monte Carlo simulations in the SU(2) and SU(3) lattice gauge 
theories showed that the color-Coulomb potential rises linearly with distance, and its string tension has 2 ~ 3 times 
larger value than that of the Wilson potential, which is an expected result from the Zwanziger 's inequality [Tt], Ell- 
in addition, it was shown that an as ymp totic scaling violation for the color-Coulomb string tension is weaker than 
that of the Wilson string tension (l9l |20|. 

In the deconfinement phase, the static potential of a quark- antiquark pair is screened due to the screening effect 
[HI H3, HH H3 ■ In contrast, it has been shown by numerical simulations that the color-Coulomb potential is a confining 
potential even in the deconfinement phase, that is, the color-Coulomb potential is not screened above the critical 
temperature of the confinement/deconfinement phase transition [17lll8|. Thus the color-Coulomb string tension does 
not serve as an order parameter for the confinement /deconfinement phase transition. This observation implies that the 
confinement is attributed to the instantaneous interaction in Coulomb gauge, whereas the confinement/deconfinement 
phase transition will be caused by the noninstantaneous interaction. 

In Gribov-Zwanziger confinement scenario, the singular behavior of the color-Coulomb potential in the infrared 
region is governed by the near-zero modes of the Faddeev-Popov (FP) operator. As Gribov pointed out, Coulomb 
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gauge does not fix a gauge completely and the gauge configurations are restricted to the Gribov region where the FP 
operator is positive. On the boundary of the Gribov region, so-called the Gribov horizon, the lowest eigenvalue of 
the FP operator vanishes. Because of entropy considerations, a probability distribution gets concentrated near the 
Gribov horizon The ghost propagator which is the expectation value of the inverse of the FP operator becomes 
singular in the infrared limit. The color-Coulomb potential in the color-singlet channel is given by 

V Coul (x-y) ee g 2 Tr[T a T b ] Q d 3 zG ac (x, z; A te )(-Vf)0°V, y; A tr )^ , (3) 

where Q is the Green's function of the FP operator and (•) denotes an Euclidean expectation value. T a (a = 1, 8) 
are the generators of su(3) Lie algebra. The singular behavior of the ghost propagator in the infrared region leads to 
the long-range interaction of the color-Coulomb potential which is responsible for the color confinement. 

Recently Greensite, Olejnik and Zwanziger studied the spectrum of the FP operator in Coulomb gauge using 
SU(2) lattice gauge simulation. The authors discussed the self-energy of an isolated quark and derived the necessary 
condition for the color confinement [26| . It was shown that the FP eigenvalue density of the lowest modes becomes 
denser as the lattice volume increases and the necessary condition is satisfied in the infinite volume limit. 

In this paper, we investigate the distribution of the FP eigenvalues in SU(3) lattice gauge simulations, and check 
whether the necessary condition for color confinement is satisfied or not. In Sec. II we discuss the necessary condition 
for the confinement in Coulomb gauge and introduce the definition of the FP operator on a lattice. Sec. Ill is devoted 
to show results of our numerical simulations. We discuss the confinement criterion also in the deconfinement phase. 
In Sec. IV, we give conclusions. 



II. COLOR-COULOMB SELF-ENERGY 

The Coulomb gauge Hamiltonian can be expressed as the sum of the gluonic part and the instantaneous part pd| : 
H = ~ [d 3 x ((E^ T (x, t)f + B?(x, tf) + \ I ' d 3 y f d 3 zp a (y, t)V ab (y, z; A^)p b (z, t). (4) 



2 J vv 1 v ' " 1 v ' ' ' 2 _ 

Here 2? l °' tr are the transverse components of the color electric field, 2?" the color magnetic field, p a (x,t) the color 
charge density. The kernel of the instantaneous interaction is given by 

V ab {y, z- A tr ) = j d 3 xG ac (y, x; A tr )(-Vl)£ cb (f , z; A tr ), (5) 

where A® ,tr are the transverse components of the gluon field. Q is the Green's function of the FP operator M ab = 
-diDf = -S ab d 2 - gf abc A c l ' tI d l whose expectation value {G ab {x, y; A u ')) = G{x-y)6 ab is the ghost propagator. The 
instantaneous interaction energy due to color charges originates from the longitudinal color electric field. In this study, 
we focus on the relation between the instantaneous interaction and the spectrum of the FP operator, and we do not 
discuss the noninstantaneous interaction which may be relevant to the confinement/deconfinement phase transition. 

The color-Coulomb self-energy for an isolated color charge, whose energy diverges in the infrared limit in a confining 
theory, is [26| 

E c = Tr[T a T b ]g 2 {V ab {x, x- A tr )). (6) 

The color-Coulomb self-energy is ultraviolet divergent in the continuum limit both in an Abelian and a nonAbelian 
gauge theories, and can be regularized by introducing the cutoff. The interesting point is that the infrared divergence 
may exist in a confining theory in the infinite volume. 

On a lattice, the FP operator is an 8V3 x 8V3 sparse matrix (V3 is the lattice 3- volume) and expressed in terms of 
SU(3) spatial link variables Ui as 



Jx,y 



- 2T b T a U l {x)5 ytX+l - 2T a T b U t (x - 1)5^ . (7) 

The Green's function of the FP operator can be expanded in terms of the eigenvectors </>"(a?) and the eigenvalues X n 
of the FP operator; 

g*(*,iM*) = y; *P ( yk (g) . (8) 
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From Eqs. J5), © and ©, we obtain 



Here Cx>(> 0) is the Casimir invariant for the representation D and F n the expectation values of the negative Laplacian 
in the FP cigcnmodes, 

F n = f d 3 xC Q (^)(-V 2 )C(f). (10) 



We define the normalized density of the FP eigenvalues 

pW = 8V3AA ' 

where N(X, X + AA) is the number of eigenvalues in the range [A, A + AA]. The total number of the eigenvalues is 8V3 
on a lattice and the FP eigenvalue density is normalized to 1. Then we have 

Wc D /™^«^», (12, 

where the upper limit of the integration X max corresponds to the UV lattice cutoff. In the Gribov-Zwanziger scenario 
the gauge configurations are restricted to the Gribov region, and therefore the lower limit of the integration is zero. 
If the condition 

hm > (13) 

A^O A 

is satisfied in the infinite volume limit, the color-Coulomb self-energy diverges in the infrared region. This is the 
necessary condition for the color confinement [26]. In this paper, we investigate whether the necessary condition for 
color confinement is satisfied in the quenched SU(3) lattice gauge simulation. 

The FP eigenvalue density of the near-zero modes is closely related to the infrared behavior of the color-Coulomb 
potential. From Eqs. © and ©, the color-Coulomb self-energy can be expressed as 

E c = V Coul (S- 2) = / ^Vcoui(p) = J ^bf Vboui(bl). (14) 

Here we introduce the ultraviolet cutoff A. If the condition (jT3"| is satisfied, the color-Coulomb self-energy diverges 
in the infrared region. Accordingly, the right-hand side of Eq. (|14| diverges in the infrared limit. It means that the 
color-Coulomb potential is more singular in the infrared region than the Coulomb potential V(p) ~ l/|p| 2 . Since 
the color-Coulomb potential provides an upper bound for the physical potential, the condition (|13[) is the necessary 
condition for the physical potential being a confining potential. 

In the Abelian gauge theory (or at the zero-th order in the coupling), the FP operator is the negative Laplacian. 
Thus the FP eigenfunctions are the plane waves and A = fc 2 . By counting the number of states in momentum space, 
it is easy to show that 



47T 2 

in the infinite volume limit. Obviously the necessary condition (|13p is not satisfied in this case 



p(X) = ^, F(X) = A, (15) 



III. NUMERICAL SIMULATIONS 



We calculate the FP eigenvalue density by the 5C/(3) lattice gauge simulations in quenched approximation. The 
lattice configurations are generated by the heat-bath Monte Carlo technique with the Wilson plaquette action (at 
(3 — 5.60 ~ 6.11) and the Iwasaki improved action (at j3 = 2.605), and we used 100 (60) thermalized configurations 
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on 8 4 ~ 20 4 (24 4 ) at zero temperature and 12 3 x 6 ~ 24 3 x 6 at finite temperature. In these simulations we adopt the 
iterative method to fix a gauge [27]. In the iterative gauge fixing procedure, we minimize the functional 

Fu[g) = ^ReTr (l - ^(x)Ui(x)g(x + i)) (16) 

with respect to the gauge transformation g(x) and find the local minimum of -Ft/ [fl 1 ] ■ The Hessian matrix associated 
with Fu\g] is the lattice FP operator in Eq. ([7]) and the local minima of the corresponding F's define the Gribov region 
[28| . Since the Hessian matrix is positive at a local minimum, the adopted method limits the lattice configurations 
to the Gribov region and the FP operator has no negative eigenvalues. We used the LAPACK package to extract 
the whole eigenvalues of the FP operator on the 8 4 lattice, while for larger lattice volumes, we used the ARPACK 
package to evaluate the lowest 1000 eigenvalues and corresponding eigenvectors of the FP operator because we are 
interested in the behavior of the low-lying FP eigenmodes. Since there are eight trivial zero modes corresponding to 
the spatially constant eigenvectors, we obtain the FP eigenvalue density from the remaining 992 eigenvalues. 

A. Full spectrum of the FP operator 



8 4 lattice, 100 confs. 



* P=6.00 




FIG. 1: The FP eigenvalue density at j3 = 6.00 on the 8 4 lattice is shown by open diamond symbols with error bars as a 
function of A. The vertical bars represent the FP eigenvalue density in the free field on the 8 4 lattice and its scaling is displayed 
in the right-hand side of the figure. The dashed line is obtained by folding the vertical bars with normalized Gauss functions 
of a fixed width to smooth out the eigenvalue distribution for comparison with the interacting case. 



The full spectrum of the FP operator on the 8 4 lattice at j3 — 6.00 is displayed in Fig. Q]in lattice units, and the 
case of the free field, Ui — /, is also shown for comparison. We show also a smooth eigenvalue distribution for the free 
field case by folding the vertical bars with normalized Gauss functions of a fixed width to smooth out the eigenvalue 
distribution for comparison with the interacting case (dashed curve). We see that the whole eigenvalues of the FP 
operator shifts to lower values; namely, the number of the lowest eigenmodes is enhanced if the interaction turns on. 

Furthermore, from this figure, we can see some bump structures, which correspond to the peaks in the case of the 
free field. The FP operator for the free field is the negative Laplacian and the FP eigenvalues are, on a L 4 lattice, 

3 

A = 4^sin 2 ^, m = 0,l,...,L-l. (17) 

i=l 

Thus, the FP eigenvalues degenerate and the eigenvalue density of the FP operator is the sum of the delta functions. The 
degeneracy of the FP eigenvalues is lost and the peaks are broadened if the interaction turns on. Eventually, the 
neighboring peaks overlap and the FP eigenvalues are distributed like those in the figure. Therefore, the appearance 
of the bumps is due to the finite volume effect. 
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FIG. 2: The FP eigenvalue density at /3 = 5.60 ~ 6.00 on the 8 4 lattice in physical units. The lattice couplings (3 = 5.60, 5.80, 6.00 
correspond to the lattice spacings a ~ 0.250, 0.150, 0.104 fm, respectively [29j . The FP eigenvalues span a wider range as the 



lattice spacing decreases since A,, 



l/o 



The FP eigenvalue density in physical units is shown in Fig. [5] The three curves deviate from each other significantly 
above A ~ 1 [GeV 2 ]. It means that the p(X) depends on the lattice cutoff seriously above A ~ 1 [GeV 2 ]. However, we 
are interested in the behavior of the FP eigenvalue density near A = and we will not discuss the cutoff dependence 
of the results anymore in this study. 



B. FP eigenvalue density with the improved action 
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FIG. 3: FP eigenvalue density on the 20 4 lattice with the Wilson action and the Iwasaki action. The lattice coupling j3 = 2.605 
for Iwasaki action corresponds to the lattice spacing a ~ 0.107 fm, and j3 = 6.00 for the Wilson action yields a ~ 0.103 fm. 



In Fig. [3] we plot the FP eigenvalue density on the 20 4 lattice with the Iwasaki action [30]. It is clear from this 
figure that the behavior of p(X) at small A with the Iwasaki action is similar to that of the Wilson action and there 
is no serious dependence of the FP eigenvalue density on the form of the lattice action. 
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C. {p(X)F(X))/X in the confined phase 
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FIG. 4: The FP eigenvalue density p(A) in the confinement phase on a variety of lattice sizes. 
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FIG. 5: The average Laplacian F(X) in the confinement phase. 

Figures 2] and [5] show p(X) and F(A) at (3 — 6.0 on 12 4 ~ 24 4 lattice volumes. We see that p(X) bends sharply as the 
lattice volume increases. On the other hand, p(X) is almost saturated above A ~ 0.15. F(X) becomes flat at smaller 
value of A as the lattice volume increases. From these figures, it seems that as A — > the FP eigenvalue density p(X) 
and the average Laplacian F(X) approach positive constants in the infinite volume limit; namely, the confinement 
criterion is satisfied in SU(3) Yang-Mills theory. 

To compare the behavior of p(X) in the nonAbelian theory with that in the Abelian theory, we plot p(X)/^/~X in 
Fig. [Sj In the Abelian theory it is constant because p(A) ~ VA (see Eq. ([To]) ). In the nonAbelian theory, we 
observe that p(X)/y/~X is almost constant above A ~ 1 [GeV 2 ]. It means that the qualitatative behaivor of p(X) in 
the nonAbelian theory is very similar to that of the Abelian theory at large A. In contrast, at small A, p{X)/yX 
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FIG. 6: p(X)/VX at (3 = 6.0 on 12 4 ~ 24 4 lattice volumes. 



is not constant and it seems to diverge as A — > in the infinite volume limit. The FP eigenvalue density in the 
nonAbelian theory shows a completely different behavior at small A compared to that of the Abelian theory and we 
see the enhancement of the near-zero modes of the FP operator. In the Gribov-Zwanziger scenario, these near-zero 
modes cause the color-Coulomb potential to be more singular in the infrared region than the simple pole. 

In the work given by Greensite et al., similar behaviors of p(X) and F(X) were observed in the SU(2) lattice gauge 
simulation. However, they did not exclude the possibility that p(X) and F(X) vanish as A — ► and they gave the 
estimates 

p(A) ~ A ' 25 , F(X) ~ A - 38 (18) 

by a scaling analysis. Then they concluded that the necessary condition for the confinement is indeed satisfied in 
SU(2) Yang-Mills theory. We have fitted the data on 24 4 lattice with the functions, 

p(X) = ciA p + c 2 Vx, F(X) = c 3 X q + c 4 A, (19) 

where we introduced the last terms c^a/A and C4A which dominate the perturbative behavior at large A. The three- 
parameter fitting gives the exponents 

p = 0.15(10), g = 0.29(4), (20) 

with x 2 /ndf = 1.60 for p(X) and x 2 /W/ = 0.763 for F(X) respectively. 

In Fig. [7] we plot p{X)F(X)/X as a function of A. As A approaches to 0, p(X)F{X)/X decreases for the free field (see 
Eq. (fTS"|) ) while increases for the interacting field. In addition, the 24 4 lattice simulation shows that p(X)F(X)/X at 
smaller A becomes flatter than the case of the 16 4 lattice simulation. From this figure, we expect that p(X)F(X)/X 
diverges or goes to positive constant, and it is unlikely that p(X)F(X)/X goes to zero as A — * in the infinite volume 
limit. Therefore, we conclude that the color-Coulomb self-energy of an isolated color charge is infrared divergent in 
SU(3) gauge theory. 

D. (p(X)F(X)) / X in the deconfined phase 

p(A), F(X) and p(X)F(X)/X in the deconfinement phase on 12 3 x6~ 24 3 x 6 lattice volumes are displayed in Figs. 
151 [91 and [T0l The simulations are carried out at (3 — 6.11 which corresponds to T/T c ~ 1.5. Here T c is the critical 
temperature of the confinement /deconfinement phase transition. These figures show that there are no drastic changes 
of these behaviors in the deconfinement phase. This is consistent with the facts that the color-Coulomb potential 
rises linearly with distance even in the deconfinement phase and the color-Coulomb string tension does not serve as 
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FIG. 7: p(A)F(A)/A vs. A in the confinement phase. 
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FIG. 8: The FP eigenvalue density p(A) in the deconfinement phase. 



an order parameter for the confinement/deconfinement phase transition. We conclude that the necessary condition 
for color confinement is satisfied even in the deconfinement phase. 

This result may not be a surprising result, since in Coulomb gauge the FP operator is a purely spatial quantity 
and does not depend on time explicitly. On the other hand, the time extent of a lattice determines temperature. 
Therefore, the FP operator is insensitive to temperature and it is natural that the spectrum of the FP operator does 
not show a critical behavior. 



IV. CONCLUSIONS 



We have calculated the eigenvalue distribution of the FP operator in Coulomb gauge using quenched SU(3) lattice 
gauge simulations. In the confinement phase, we observe the accumulation of the near-zero eigenvalues of the FP 
operator at large lattice volumes. Moreover, the lattice simulations reveal that the average Laplacian approaches 
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FIG. 9: The average Laplacian F(X) in the decomfinement phase. 
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FIG. 10: p(X)F(X)/X vs. A in the deconfinement phase. 



positive constant as A — > 0. We conclude that the confinement criterion is satisfied in the SU(3) gauge theory. This 
supports the Gribov-Zwanziger confinement scenario. The results we obtained are qualitatively consistent with those 
of the SU(2) lattice simulation carried out by Greensite et al. 

The near-zero modes of the FP operator survive above the critical temperature, and the behaviors of the FP eigen- 
value density and the average Laplacian in the deconfinement phase are qualitatively the same as in the confinement 
phase. Accordingly, the confinement criterion is satisfied even in the deconfinement phase in SU(3) gauge theory. It 
is not surprising that the spectrum of the FP operator is insensitive to temperature, because the FP operator is a 
spatial quantity. We note that the criterion is not a sufficient condition but a necessary condition for the confinement; 
namely, the color-Coulomb energy is not the ground state energy but the excited state energy of color charges [l4l |. 
If we take the noninstantaneous interaction into account when discussing in the deconfinement phase, the energy of 
an isolated color charge will be finite in the infrared limit and the isolated color charge can exist. 

As we have shown in this paper, the spectrum of the FP operator does not change drastically above the critical 
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temperature. This would indicate that confining features survive even in the deconfincment phase. Actually, it is 
known that the spatial Wilson loop which is a gauge invariant quantity shows the area law behavior above the critical 
temperature. Therefore, we expect that further studies in Coulomb gauge provide insight into the understanding of 
the strongly correlated quark-gluon plasma. 

The color-Coulomb potential can be obtained by calculating the FP eigenvalues and eigenvectors. It is valuable 
to see whether the lowest eigenmodes of the FP operator produce the linearly rising behavior of the color-Coulomb 
potential for large quark separations. We address this issue in our future investigation. 
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